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1.  NTRODUCTION 


In  a  collaboration  with  Dr.  Tonia  Rex  we  developed  a  mouse  model  of  blast  injury 
to  the  eye,  which  accurately  mimics  the  traumatic  blast  injury  increasingly 
suffered  by  warriors  under  current  battlefield  conditions  (Hines-Beard  et  al., 
2012).  Using  this  mouse  model  in  combination  with  a  powerful  combination  of 
systems  biology,  microarray  analysis,  expression  genetics,  and  bioinformatics, 
we  are  defining  the  genetic  networks  activated  by  the  ocular  blast  injury.  At  the 
heart  of  our  approach  is  a  genetic  reference  panel  of  mice,  the  unique  resource 
of  BXD  recombinant  inbred  (Rl)  strain  set.  The  set  of  Rl  strains  was  produced 
from  a  genetic  cross  between  the  C57BU6J  mouse  and  the  DBA/2J  mouse. 
Using  60  BXD  strains  provides  a  new  and  powerful  method  to  defining  elements 
in  the  genome  regulating  the  response  of  the  eye  to  blast  injury.  This  allows  us  to 
generate  specific,  testable  hypotheses  to  define  the  pathways  that  regulate  the 
response  of  the  eye  to  blast  injury  and  reactive  responses  in  the  retina.  As  more 
diverse  gene  expression  data  sets  become  available,  comparison  of  gene 
expression  and  regulation  in  different  biological  contexts  will  help  identify  the 
regulatory  elements  controlling  the  injury  response  of  the  eye  and  the  retina.  We 
are  currently  identifying  genetic  networks  activated  by  blast  injury  and  the 
genomic  regions  controling  these  networks.  We  are  also  identifying  markers  for 
retinal  injury  and  potential  targets  for  intervention. 


2  KEYWORDS 
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3  ACCOMPLISHMENTS 


Maior  Goals: 

Task  1)  Quantify  the  strain-to-strain  differences  in  the  severity  of  blast- 
induced  ocular  pathologies,  using  a  set  of  60  BXD  Rl  mouse  strains  and 
map  the  genomic  loci  that  regulate  the  response  of  the  eye  to  blast  injury. 
In  this  Task  we  were  measuring  intraocular  pressure  (IOP),  central  corneal 
thickness  (CCT)  and  visual  acuity. 

Task  2)  Define  the  genetic  networks  activated  by  blast  injury  in  the  eye  and 
in  the  retina,  using  transcriptome-wide  profiling  across  the  BXD  Rl  strain 
set.  We  are  using  the  Affymetrix  GeneChip  Gene  2.0  ST  Mouse  Array  to 
characterize  the  changes  occurring  following  a  blast  injury  to  the  eye  in  60  BXD 
strains.  There  were  several  major  benefits  to  using  the  new  Affymetrix  array. 
Specifically,  there  are  probes  for  7,000  non-coding  RNAs  (RNA  that  is  not 
converted  to  protein  but  does  affect  the  functioning  of  the  cell).  We  are  now 
finding  out  that  many  of  these  non-coding  RNAs  play  extremely  important  roles  in 
the  body.  Within  these  7,000  probes,  588  encode  microRNAs  (small  RNAs  that 
regulate  protein  expression).  We  are  creating  an  entire  normal  retina  dataset 
using  the  Affymetrix  GeneChip  Gene  2.0  Mouse  Array  and  comparing  this  data 
set  to  a  dataset  from  retinas  5  days  after  a  50psi  blast  injury  to  the  eye. 

Task  3)  Define  biomarkers  that  can  predict  the  severity  of  injury  and 
eventual  outcomes. 

This  portion  of  our  study  was  to  begin  in  the  latter  years  of  the  grant  (Months  40 
to  48).  We  are  using  this  to  characterize  the  50-psi  blast  injury  in  advance  of 
resuming  the  blast  microarray  study  on  the  BXD  Rl  strain  set.  Immunostaining 
sections  of  retina  revealed  that  SOX  11  was  upregulated  in  the  neurons  of  the 
inner  retina  following  blast.  SOX  11  labeled  cells  in  the  ganglion  cell  layer  and  the 
inner  nuclear  layer.  In  the  ganglion  cell  layer  SOX1 1  labeled  a  majority  of  the 
cells,  indicating  that  it  was  labeling  most  ganglion  cells  and  displaced  amacrine 
cells.  Once  the  datasets  are  fully  implemented,  we  will  be  able  to  accurately 
define  the  changes  occurring  within  the  injured  retina. 

Accomplishments  Under  These  Goals: 

Taskl: 

At  the  present  time  we  have  measured  IOP  and  central  corneal  thickness  on  33 
strains  of  mice  before  and  after  a  50psi  blast  injury  to  the  eye.  At  the  present  time 
there  are  117  animals  total  in  the  dataset.  When  we  run  a  student  t  test  on  the 
data  there  was  no  significant  difference  in  CCT  or  IOP  before  and  after  blast  in 
the  control  eyes.  This  is  expected.  We  also  did  not  see  a  significant  difference  in 
either  CCT  or  IOP  5  days  after  a  50psi  blast  to  the  experimental  eye.  This  is 
unexpected.  We  intend  to  continue  measuring  CCT  and  IOP  measurements  for 
the  next  quarter. 
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Task  2. 


A)  We  have  completed  the  construction  of  the  DoD  CDMRP  Normal  Retinal 
Dataset.  Using  the  Affymetrix  Mouse  Gene  2.0  ST  array,  to  interrogate  all  exons 
of  traditional  protein  coding  genes,  non-coding  RNAs  and  microRNAs.  These 
data  are  presented  in  a  highly  interactive  database  within  the  GeneNetwork 
website.  In  the  Normal  Retina  Database,  we  quantified  mRNA  levels  of  the 
transqriptome  from  retinas  using  the  Affymetrix  Mouse  Gene  2.0  ST  array.  The 
Normal  Retina  Database  consists  of  gene  expression  data  from  male  and  female 
mice.  The  dataset  includes  a  total  of  55  BXD  Rl  strains,  the  parental  strains 
(C57BI/6J  and  DBA/2J),  and  a  reciprocal  cross.  In  combination  with 
GeneNetwork,  the  DoD  (Department  of  Defense)  CDMRP  (Congressionally 
Directed  Medical  Research  Programs)  Normal  Retina  Database  provides  a  large 
resource  for  mapping,  graphing,  analyzing,  and  testing  complex  genetic  networks. 
Protein-coding  and  non-coding  RNAs  can  be  used  to  map  quantitative  trait  loci 
(QTLs)  that  contribute  to  expression  differences  among  the  BXD  strains  and  to 
establish  links  between  classical  ocular  phenotypes  associated  with  differences 
in  genomic  sequence.  With  this  resource  we  are  able  to  extract  transcriptome 
signatures  for  retinal  cells  and  to  define  genetic  networks  associated  with  the 
maintenance  of  the  normal  retina.  Ultimately,  we  will  use  this  database  to  define 
changes  occurring  following  blast  injury  to  the  retina.  The  DoD  CDMRP  Normal 
Retina  Database  uses  the  Affymetrix  MouseGene  2.0  ST  Array  (May  15  2015). 
The  RMA  analysis  and  scaling  was  conducted  by  Arthur  Centeno.  This  data  set 
consists  of  55  BXD  strains,  C57BU6J,  DBA/2J,  an  FI  cross  between  C57BU6J 
and  DBA/2J.  A  total  of  58  strains  were  quantified.  There  is  a  total  of  222 
microarrays.  All  of  the  data  from  each  of  the  microarrays  used  in  this  dataset  is 
publically  available  on  GeneNetwork.org. 

Mice  were  killed  by  rapid  cervical  dislocation.  Retinas  were  removed  immediately 
and  placed  in  1  ml  of  160  U/ml  Ribolockfor  1  min  at  room  temperature.  The 
retinas  were  removed  from  the  eye  and  placed  in  Hank's  Balanced  Salt  solution 
with  RiboLock  in  50p1  Ribolock  (Thermo  Scientific  RiboLock  RNase#E00381 
40U/pl2500U)  and  stored  in-80°C.  The  RNAwas  isolated  using  aQiaCube.  All 
RNA  samples  were  checked  for  quality  before  running  microarrays.  The  samples 
were  analyzed  using  the  Agilent  2100  Bioanalyzer.  The  RNA  integrity  values 
ranged  from  7.0  to  1 0.  Our  goal  was  to  obtain  data  for  independent  biological 
sample  pools  from  both  sexes  for  most  lines  of  mice.  The  four  batches  of  arrays 
included  inthis  final  data  set  collectively  represent  a  reasonably  well-balanced 
sample  of  males  and  females,  in  general  without  within-strain-by-sex  replication. 
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The  data  is  presented  using  the  Affymetrix  Mouse  Gene  2.0  ST  Array.  These 
expression  arrays  have  been  designed  with  a  median  of  22  unique  probes  per 
transcript.  Each  unique  probe  is  25  bases  in  length,  which  means  that  the  array 
measures  a  median  of  550  bases  per  transcript.  The  arrays  provide 
comprehensive  transcriptome  coverage  with  over  30,000  coding  and  non-coding 
transcripts.  In  addition  there  is  coverage  for  over  600  microRNAs.  For  some 
arrays  the  RNA  was  pooled  from  two  retinas  and  other  arrays  were  run  on  a 
single  retina.  Dr.  XiangDi  Wang  (UTHSC)  and  Rebecca  King  (Emory)  were 
involved  in  the  retinal  extractions  and  isolation  of  RNA.  Two  different  research 
cores  ran  the  Affymetrix  arrays:  the  Molecular  Resource  Center  at  UTHSC  (Dr. 
William  Taylor  Director)  and  the  Integrated  Genomics  Core  at  Emory  University 
by  Robert  B.  Isett  (Dr.  Michael  E.  Zwick,  Director).  In  a  separate  set  of 
experiments  we  tested  a  set  of  arrays  from  C57BU6J  retinas  run  at  each  facility 
to  determine  if  there  were  batch  effects  or  other  confounding  differences  inthe 
results.  We  could  not  detect  any  significant  difference  inthe  arrays  run  at 
UTHSC  or  at  Emory  University.  Thus,  we  have  included  both  sets  of  data  into  the 
analysis. 

Publication:  Rebecca  King,  Lu  Lu,  Robert  W.  Williams  and  Eldon  E.  Geisert, 
(2015)  Gene  Expression  and  Genetic  Networks  inthe  Mouse  Retina,  Submitted 
to  Molecular  Vision  (under  review). 


Task  3)  We  have  identified  a  list  of  potential  biomarkers  for  injury  to  the  retinal 
ganglion  cells.  The  best  marker  is  SOX1 1  (manuscript  being  revised).  We  are 
using  this  to  characterize  the  50psi  blast  injury  in  advance  of  resuming  the  blast 
microarray  study  on  the  BXD  Rl  strain  set.  Immunostaining  sections  of  retina 
revealed  that  SOX  11  was  upregulated  inthe  neurons  of  the  inner  retina  following 
blast.  SOX1 1  labeled  cells  in  the  ganglion  cell  layer  and  the  inner  nuclear  layer. 
In  the  ganglion  cell  layer  SOX1 1  labeled  a  majority  of  the  cells,  indicating  that  it 
was  labeling  most  ganglion  cells  and  displaced  amacrine  cells.  Amacri  ne  cells  in 
the  inner  nuclear  layer  were  also  lightly  labeled  by  SOX11.  On  immunoblots 
there  was  approximately  a  2-fold  increase  in  the  intensity  of  the  SOX11  band. 
The  manuscript  describing  these  results  is  currently  being  revised  to  respond  to 
the  editorial  suggestions. 


Publication:  Felix  L.  Struebing,  Steven  G.  Hart,  Joseph  M.  Caron,  XiangDi  Wang 
and  Eldon  E.  Geisert  (2015)  Upregulation  of  SOX  11  in  Retinal  Ganglion  Cells 
Following  Injury.  Molecular  Vision  (under  revision  priorto  resubmission). 


Training  and  Professional  Development  Opportunities: 
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Nothing  to  Report 
Dissemination  of  Results: 

Invited  Talks: 

2014  Genetic  Network  of  Innate  Immunity  in  the  Retina:  Relevance  to  CNS  Injury 
and  Disease  Department  of  Molecular  Physiology  and  Biophysics  University  of 
Iowa. 

2014  Genetic  Network  of  Innate  Immunity  in  the  Retina:  Relevance  to  CNS  Injury 
and  Disease  VA  Atlanta  GA. 


ARVO  Meeting: 

Sidhu  C,  Lyuboslavsky  P,  Chrenek  MA,  Struebing  FL,  Sellers  JT,  Setterholm  NA, 
McDonals  FE,  Boatright  JH,  Geisert  EE,  luvone  PM:  Traumatic  Blast-Induced 
Closed  Globe  Injury  Reduces  Visual  Function  in  Retinal  Ganglion  Cells  of  Thyl- 
CFP  mice:  Mitigation  by  a  Small  Molecule  TrkB  Activator",  The  Association  for 
Research  in  Vision  and  Ophthalmology  (ARVO)  annual  meeting,  Denver  201 5. 


Papers  Submitted  for  Publication: 

Rebecca  King,  Lu  Lu,  Robert  W.  Williams  and  Eldon  E.  Geisert,  (2015)  Gene 
Expression  and  Genetic  Networks  in  the  Mouse  Retina,  Submitted  to  Molecular 
Vision  (under  review). 

Felix  L.  Struebing,  Steven  G.  Hart,  Joseph  M.  Caron,  XiangDi  Wang  and  Eldon  E. 
Geisert  (2015)  Upregulation  of  SOX1 1  in  Retinal  Ganglion  Cells  Following  Injury. 
Molecular  Vision  (under  revision  prior  to  resubmission). 


Plans  for  Next  Reporting  Period  to  Accomplish  the  Goals: 

1)  Finish  Beta  testing  and  open  the  DoD  CDMRP  Normal  Retina  Database  to  the 
public  in  August  2015. 

2)  We  will  complete  the  DoD  Blast  Microarray  Database.  This  will  entail  running 
several  hundred  microarrays  along  with  the  costs  of  analysis.  We  anticipate  this 
to  consume  a  significant  amount  of  the  remaining  budget. 

3)  Prepare  manuscript  describing  the  Blast  database. 

4)  Release  the  Blast  database  to  the  public. 
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4.  IMPACT 


Impact  on  the  Development  of  the  Principal  Discipline  of  the  Project: 

Once  the  proposed  studies  are  completed  they  will  provide  a  comprehensive 
analysis  of  the  molecular  pathways  activated  in  the  retina  by  blast  injury  to  the 
eye. 


Impact  on  Other  Disciplines: 

When  developing  Biomarkers  for  retinal  injury,  our  microarray  dataset  will  provide 
a  means  to  determine  if  any  specific  biomarker  could  have  originated  from  the 
retinal  injury  itself. 

Impact  on  Society  Beyond  Science  and  Technology: 

Nothing  to  Report 
5)  Changes/Problems 

Changes  in  Approach  and  Reasons  for  Change: 

None 

Actual  or  Anticipated  Problems  or  Delays  and  Actions  or  Plans  to  Resolve 
Them: 

None 

Changes  that  had  Significant  Impact  on  Expenditures: 

None 


Significant  Changes  in  the  Use  or  Care  of  Human  Subjects  Vertebrate 
Animals  Biohazards,  or  Select  Agents: 

None 


6.  PRODUCTS 

Publications,  conference  papers,  and  presentations: 

2014  Genetic  Network  of  Innate  Immunity  in  the  Retina:  Relevance  to  CNS  Injury 
and  Disease  Department  of  Molecular  Physiology  and  Biophysics  University  of 
Iowa. 
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201 4  Genetic  Network  of  Innate  Immunity  in  the  Retina:  Relevance  to  CNS  Injury 
and  Disease  VA  Atlanta  GA. 


ARVO  Meeting: 

Sidhu  C,  Lyuboslavsky  P,  Chrenek  MA,  Struebing  FL,  Sellers  JT,  Setterholm  NA, 
McDonals  FE,  Boatright  JH,  Geisert  EE,  luvone  PM:  "Traumatic  Blast-Induced 
Closed  Globe  Injury  Reduces  Visual  Function  in  Retinal  Ganglion  Cells  of  Thyl- 
CFP  mice:  Mitigation  by  a  Small  Molecule  TrkB  Activator",  The  Association  for 
Research  in  Vision  and  Ophthalmology  (ARVO)  annual  meeting,  Denver  201 5. 


Papers  Submitted  for  Publication: 

Rebecca  King,  Lu  Lu,  Robert  W.  Williams  and  Eldon  E.  Geisert,  (2015)  Gene 
Expression  and  Genetic  Networks  in  the  Mouse  Retina,  Submitted  to  Molecular 
Vision  (under  review). 

Felix  L.  Struebing,  Steven  G.  Hart,  Joseph  M.  Caron,  XiangDi  Wang  and  Eldon  E. 
Geisert  (2015)  Upregulation  of  SOX1 1  in  Retinal  Ganglion  Cells  Following  Injury. 
Molecular  Vision  (under  revision  prior  to  resubmission). 

Website(s)  or  other  Internet  site(s): 

The  DoD  CDMRP  Retina  Affy  MoGene  2.0  ST  Database  and  the  DoD  TATRC 
Retina  Affy  MoGene  2.0  ST  Exon  Level  Database  are  hosted  on 
GeneNetwork.org.  This  database  will  be  open  to  the  public  in  August  2015. 
These  datasets  describe  gene  expression  in  the  normal  retina  of  the  BXD  Strains. 
Both  databases  can  be  found  under  Mice,  BXD,  retina  and  then  either  DoD 
CDMRP  Retina  Affy  MoGene  2.0  ST  Database  or  DoD  CDMRP  Retina  Affy 
MoGene  2.0  ST  Exon  Level  Database. 


Technologies  or  techniques: 

None 

Inventions,  patent  applications,  and/or  licenses: 

None 

Other  products: 

None 

7. PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
What  individuals  have  worked  on  the  project? 
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At  Emory  University  (7/1 5/1 4 to  present): 

Becky  King,  Research  Technician  (50%  effort) 

Eldon  E.  Geisert,  Principal  Investigator  (25%  effort) 

We  have  been  collaborating  with  Dr.  Mike  luvone  to  construct  and  test  a  new 
blast  gun.  We  are  currently  in  the  process  of  writing  a  manuscript  describing  the 
effects  of  a  50psi  blast  to  the  mouse  eye. 

Has  there  been  a  change  in  the  other  active  support  of  the  PD/PI(s)  or 
senior/key  personnel  since  the  last  reporting  period? 

No. 

What  other  organizations  have  been  involved  as  partners? 

None 

8.  SPECIAL  REPORTING  REQUIREMENTS 
None 


9.APPENDI  CES 

A)  Preprint  of  the  DoD  CDMRP  Normal  Retina  Database  paper. 


8)  Preprint  of  the  SOX1 1  (marker  of  retinal  injury)_  paper. 
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30  Abstract 

31 

32  Purpose:  Differences  in  gene  expression  provide  diverse  retina  phenotypes  and 

33  may  also  contribute  to  susceptibility  to  injury  and  disease.  The  present  study 

34  defines  the  transcriptome  of  the  retina  in  the  BXD  Rl  strain  set.  Using  the 

35  Affymetrix  Mouse  Gene  2.0  ST  array,  to  interrogate  all  exons  of  traditional  protein 

36  coding  genes,  non-coding  RNAs  and  microRNAs.  These  data  are  presented  in  a 

37  highly  interactive  database  within  the  GeneNetwork  website. 

38 

39  Methods:  In  the  Normal  Retina  Database,  quantified  mRNA  levels  of  the 

40  transcriptome  from  retinas  using  the  Affymetrix  Mouse  Gene  2.0  ST  array.  The 

41  Normal  Retina  Database  consists  of  gene  expression  data  from  male  and  female 

42  mice.  The  dataset  includes  a  total  of  55  BXD  Rl  strains,  the  parental  strains 

43  (C57BI/6J  and  DBA/2J),  and  a  reciprocal  cross. 

44 

45  Results:  In  combination  with  GeneNetwork,  the  DoD  (Department  of  Defense) 


46  CDMRP  (Congressionally  Directed  Medical  Research  Programs)  Normal  Retina 

47  Database  provides  a  large  resource  for  mapping,  graphing,  analyzing,  and 

48  testing  complex  genetic  networks.  Protein-coding  and  non-coding  RNAs  can  be 

49  used  to  map  quantitative  trait  loci  (QTLs)  that  contribute  to  expression 

50  differences  among  the  BXD  strains  and  to  establish  links  between  classical 

51  ocular  phenotypes  associated  with  differences  in  genomic  sequence.  With  this 

52  resource  we  are  able  to  extract  transcriptome  signatures  for  retinal  cells  and  to 

53  define  genetic  networks  associated  with  the  maintenance  of  the  normal  retina. 

54  Ultimately,  we  will  use  this  database  to  define  changes  occurring  following  blast 

55  injury  to  the  retina. 

56 

57  Conclusions:  The  high  level  of  variation  in  mRNA  levels  found  among  BXD  Rl 

58  strains  of  mice  make  it  possible  to  identify  expression  networks  underlying 

59  differences  in  retina  structure  and  function. 

60 
61 
62 

63 

64 

65 

66 

67  INTRODUCTION 

68  Large-scale  sequencing  initiatives  have  led  to  a  new  era  in  understanding 

69  gene  and  genome  functions  [1  -5].  There  is  now  an  acute  need  for  powerful 

70  approaches  to  integrate  and  analyze  massive  omics  data  sets.  For  vision 

71  research  there  are  many  known  single  gene  variants  that  cause  vision  loss, 

72  including  retinitis  pigmentosa  [6-9],  Usher’s  syndrome  [10,  11]  and  some  forms  of 

73  glaucoma  [1 2],  However,  many  ocular  diseases  have  a  complex  genetic  basis 

74  with  multiple  chromosomal  loci  contributing  to  differences  in  susceptibility  and 

75  severity  of  disease.  Two  prominent  examples  are  glaucoma  [1 3-1 5]  and  age- 

76  related  macular  degeneration  [16,  17].  In  addition,  the  response  of  the  eye  and 

77  retina  to  trauma  is  driven  by  a  host  of  different  genes  expressed  in  a  large 

78  number  of  different  cell  types. 

79 

80  Until  recently,  it  was  extremely  difficult  to  define  the  genetic  and  molecular  basis 

81  of  complex  diseases  or  to  adequately  monitor  the  response  of  eye  and  retina  to 

82  injury.  We  have  usee  a  novel  and  powerful  approach  relying  on  systems  biology 

83  and  a  mouse  genetic  reference  panel,  the  BXD  family  of  recombinant  inbred  (Rl) 

84  strains.  This  resource  is  particularly  well  suited  to  define  complex  genetic 

85  networks  that  are  also  active  in  human  diseases.  This  approach  not  only  allows 

86  us  to  identify  specific  gene  variants  involved  in  the  retinal  disease  and  response 

87  to  injury,  but  also  allows  us  to  place  corresponding  moelcular  changes  in  a  global 

88  contest  in  eye  and  retina. 

89 

90 


91  Our  initial  efforts  of  our  group  explored  the  genetic  diversity  of  the  BXD 

92  family  of  strains  to  define  genetic  networks  active  in  the  eye  (see  data  sets  and 

93  refs  [18])  and  ([19]). 

94  In  this  study  we  have  created  the  new  mouse  retinal  database  that  offers  a 

95  more  complete  description  of  the  mouse  transcriptome.  This  resource  uses  the 

96  genetic  covariance  of  expression  across  a  panel  of  55  BXD  strains  to  identify 

97  cellular  signatures  and  genetic  networks  within  the  mouse  retina.  The  array  we 

98  have  used  provides  expression  profiling  at  the  exon  level  for  26,1 91  well- 

99  established  annotated  transcripts,  as  well  as  9,049  non-coding  RNAs  including 

100  over  600  microRNAs.  Using  the  bioinformatics  tools  located  within  GeneNetwork 

101  (genenetwork.org),  we  examine  a  cellular  signature  of  RPE  cells.  We  also  have 

102  analyzed  a  genetic  and  molecular  network  involved  in  neuronal  development  and 

103  axon  growth.  In  both  of  these  examples  we  highlight  the  specific  benefits  of  the 

104  new  database  with  a  special  emphasis  on  microRNAs,  non-coding  RNAs  and  the 

105  exon  level  data  available  with  the  Affymetrix  MouseGene  2.0  ST  array. 

106 

107  MATERIALS  and  METHODS 

108 

109  The  DoD  (Department  of  Defense)  CDMRP  (Congressionally  Directed  Medical 

110  Research  Programs)  Normal  Retina  Database  uses  the  Affymetrix  MouseGene 

111  2.0  ST  Array  (May  15  2015).  The  RMA  analysis  and  scaling  was  conducted  by 

112  Arthur  Centeno.  This  data  set  consists  of  55  BXD  strains,  C57BL/6J,  DBA/2J,  an 

113  FI  cross  between  C57BL/6J  and  DBA/2J.  A  total  of  58  strains  were  quantified. 

114  There  is  a  total  of  222  microarrays.  All  of  the  data  from  each  of  the  microarrays 

115  used  in  this  dataset  is  publically  available  on  GeneNetwork.org. 

116 

117  This  is  RMA  expression  data  that  has  been  normalized  using  what  we  call  a  2z+8 

118  scale,  but  without  special  correction  for  batch  effects.  The  data  for  each  strain 

119  was  computed  as  the  mean  of  four  samples  per  strain.  Expression  values  on  a 

120  log2  scale  range  from  3.81  to  14.25  (10.26  units),  a  nominal  range  of 

121  approximately  1 ,000-fold.  After  taking  the  log2  of  the  original  non-logged 

122  expression  estimates,  we  convert  data  within  an  array  to  a  z-score.  We  then 

123  multiply  the  z-score  by  2.  Finally,  we  add  8  units  to  ensure  that  no  values  are 

124  negative.  The  result  is  a  scale  with  a  mean  expression  of  the  probes  on  the  array 

125  of  8  units  and  a  standard  deviation  of  2  units.  A  two-fold  difference  in  expression 

126  is  equivalent  roughly  to  1  unit  on  this  scale.  The  lowest  level  of  expression  is  3.81 

127  ( Olfrl  186)  from  DoD  CDMRP  (Normal  Retina  Database  uses  the  Affymetrix 

128  MouseGene  2.0  ST  Array  (May  15  2015)  The  highest  level  of  expression  is 

129  Rhodopsin  for  17462036  (Rho).  Highest  single  value  is  about  14.25. 

130 

131 

132  About  the  cases  used  to  generate  this  set  of  data: 

133  Almost  all  animals  are  young  adults  between  60  and  100  days  of  age.  We 

134  measured  expression  in  conventional  inbred  strains,  BXD  recombinant  inbred 

135  (Rl)  strains,  and  reciprocal  FIs  between  C57BL/6J  and  DBA/2J. 

136 


137  BXD  strains: 

138  The  first  32  of  these  strains  are  from  the  Taylor  series  of  BXD  strains  generated 

139  at  the  Jackson  Laboratory  by  Benjamin  A.  Taylor.  BXD1  through  BXD32  were 

140  started  in  the  late  1970s,  whereas  BXD33  through  42  were  started  in  the  1990s. 

141  BXD43  and  higher  were  bred  by  Lu  Lu,  Jeremy  Peirce,  Lee  M.  Silver,  and  Robert 

142  W.  Williams  starting  in  1997  using  B6D2  generation  10  advanced  intercross 

143  progeny.  This  modified  breeding  protocol  doubles  the  number  of  recombinations 

144  per  BXD  strain  and  improves  mapping  resolution  [20].  All  of  the  Taylor  series  of 

145  BXD  strains  and  many  of  the  new  BXD  strains  are  available  from  the  Jackson 

146  Laboratory.  Several  strains  were  specifically  excluded  from  the  dataset.  For  the 

147  BXD43  and  higher,  the  DBA/2J  parent  carried  both  the  Tyrp-1  mutation  and  the 

148  Gpnmb  mutation  and  these  two  mutations  produce  pigment  dispersion  glaucoma. 

149  All  of  the  mice  carrying  these  two  mutations  were  not  included  in  the  dataset: 

150  BXD53,  BXD55,  BXD62,  BXD66,  BXD68,  BXD74,  BXD77,  BXD81 ,  BXD88, 

151  BXD89,  BXD95  and  BXD98.  In  addition  BXD24  was  omitted,  since  it  developed  a 

152  spontaneous  mutation,  rd16  (Cep290)  which  resulted  in  retinal  degeneration  and 

153  was  renamed  BXD24b/TyJ  [21].  Several  additional  strains  were  excluded  due  to 

154  abnormally  high  Gfap  levels  observed  in  our  Full  HEI  Retina  (April  201 0)  dataset, 

155  these  include:  BXD32,  BXD49,  BXD70,  BXD83  and  BXD89. 

156 

157  Tissue  preparation  protocol. 

158  Mice  were  killed  by  rapid  cervical  dislocation.  Retinas  were  removed  immediately 

159  and  placed  in  1  ml  of  160  U/ml  Ribolockfor  1  min  at  room  temperature.  The 

160  retinas  were  removed  from  the  eye  and  placed  in  Hank’s  Balanced  Salt  solution 

161  with  RiboLock  in  50pl  RiboLock  (Thermo  Scientific  RiboLock  RNase  #EO0381 

162  40U/pl  2500U)  and  stored  in  -80 °C.  The  RNA  was  isolated  using  a  QiaCube  and 

163  the  in  column  DNAse  procedure.  All  RNA  samples  were  checked  for  quality 

164  before  running  microarrays.  The  samples  were  analyzed  using  the  Agilent  2100 

165  Bioanalyzer.  The  RNA  integrity  values  for  ranged  from  7.0  to  1 0.  Our  goal  was  to 

166  obtain  data  for  independent  biological  sample  pools  from  both  sexes  for  most 

167  lines  of  mice.  The  four  batches  of  arrays  included  in  this  final  data  set  collectively 

168  represent  a  reasonably  well-balanced  sample  of  males  and  females,  in  general 

169  without  within-strain-by-sex  replication. 

170 

171 

172  Affymetrix  Mouse  Gene  2.0  ST  Array:  These  expression  arrays  have  been 

173  designed  with  a  median  of  22  unique  probes  per  transcript.  Each  unique  probe  is 

174  25  bases  in  length,  which  means  that  the  array  measures  a  median  of  550  bases 

175  per  transcript.  The  arrays  provide  comprehensive  transcriptome  coverage  with 

176  over  30,000  coding  and  non-coding  transcripts.  In  addition  there  is  coverage  for 

177  over  600  microRNAs.  For  some  arrays  the  RNA  was  pooled  from  two  retinas  and 

178  for  other  arrays  were  run  on  a  single  retina.  Dr.  XiangDi  Wang  (UTHSC)  and 

179  Rebecca  King  (Emory)  were  involved  in  the  retinal  extractions  and  isolation  of 

180  RNA.  The  Affymetrix  arrays  were  run  by  two  different  research  cores:  the 

181  Molecular  Resource  Center  at  UTHSC  (Dr.  William  Taylor  Director)  and  the 

182  Integrated  Genomics  Core  at  Emory  University  by  Robert  B.  Isett  (Dr.  Michael  E. 


183  Zwick,  Director).  In  a  separate  set  of  experiments  we  tested  a  set  of  arrays  from 

184  C57BL/6J  retinas  run  at  each  facility  to  determine  if  there  were  batch  effects  or 

185  other  confounding  differences  in  the  results.  We  could  not  detect  any  significant 

186  difference  in  the  arrays  run  at  UTHSC  or  at  Emory  University.  Thus,  we  have 

187  included  both  sets  of  data  into  the  analysis. 

188 

189 

190 

191  RESULTS 

192 

193 

194  The  DoD  CDMRP  Retina  Database  presents  the  retinal  transcriptome  profiles  of 

195  52  BXD  Rl  strains  in  a  highly  interactive  website,  GeneNetwork.  There  are  two 

196  separate  presentations  of  the  microarray  data.  The  first  is  at  the  gene  level  (DoD 

197  CDMRP  Retina  Affy  MoGene  2.0  ST  (May15)  RMA  Gene  Level  Database)  and 

198  the  same  data  presented  at  the  exon  level  (DoD  CDMRP  Retina  Affy  MoGene 

199  2.0  ST  (May15)  RMA  Exon  Level  Database).  For  the  analysis  of  these  dataset 

200  there  is  a  suite  of  bioinformatics  tools  integrated  into  the  GeneNetwork  website. 

201  These  tools  allow  for:  the  identification  of  genes  that  vary  across  the  BXD  Rl 

202  strains,  the  construction  of  genetic  networks  controlling  the  development  of  the 

203  mouse  retina,  and  the  identification  of  genomic  loci  underlying  complex  traits  in 

204  the  retina.  In  the  present  paper  we  present  these  two  new  datasets  and  illustrate 

205  their  use  with  two  examples.  The  first  was  to  identify  genetic  signatures  of  the 

206  retinal  pigment  epithelium  (RPE).  The  second  will  identify  a  genetic  network 

207  associated  with  roundabout  homolog  2  ( Robo2 )  gene  and  the  modulating  axonal 

208  growth. 

209 

210  Cellular  Signature  of  RPE  in  the  DoD  CDMRP  Retina  Database 

211 

212  The  DoD  CDMRP  Retina  Database  has  a  unique  signature  for  RPE  cells.  When 

213  looking  at  the  expression  of  the  RPE  marker  Rpe65  there  was  an  almost  biphasic 

214  distribution  of  expression  (Figure  1 ).  Many  of  the  strains  expressed  relatively  low 

215  levels  of  Rpe65  (approximately  7  units  on  our  scale)  while  other  strains  had  high 

216  levels  of  expression  ranging  from  2  to  8  fold  higher  (8  to  1 1  units).  When  we 

217  examined  the  dataset  for  genes  with  similar  expression  across  the  BXD  strains,  a 

218  list  of  genes  uniquely  expressed  in  RPE  was  observed  (Table  1 ).  This  cellular 

219  signature  represents  genes  that  are  uniquely  expressed  within  the  RPE, 

220  including:  Rgr  (retinal  G  protein  coupled  receptor),  Lrat  (Lecithin-retinol 

221  acyltransferase),  Rdh5  (retinol  deydrogenase  5),  Trf  (transferrin)  and  Rrh  (retinal 

222  pigment  epithelium  derived  rhodopsin  homolog).  This  signature  can  also  be 

223  thought  of  as  the  result  of  genetic  networks  that  drive  gene  expression  within  a 

224  given  cell  type.  With  the  new  Affymetrix  chip  we  not  only  have  protein-coding 

225  genes  that  correlate  with  Rpe65,  but  we  also  have  microRNAs  and  non-coding 

226  RNAs.  If  we  examine  the  top  500  correlates  of  Rpe65  (all  of  which  have  a 

227  correlation  higher  than  0.8  with  Rpe65),  there  are  five  microRNAs  present:  Mir98, 

228  Mir666,  Mir449a,  Mir301b  and  Mir28b.  Using  the  bioinformatics  tools  on 


229  TargetScan  (Targetscan.org)  [22-24]  we  were  able  to  predict  targets  for  each  of 

230  the  microRNAs  from  the  top  500  correlates  of  Rpe65.  One  microRNA,  Mir666 

231  did  not  appear  on  the  Targetscan  website.  The  remaining  4  microRNAs  did 

232  appear  on  the  website.  When  scanned  for  targets  Mir98  had  29  targets  in  the 

233  RPE  signature,  Mir449a  hadl  4  targets,  Mir301b  had  1 3  targets  and  Mir28b  had  1 

234  target.  This  type  of  analysis  may  be  one  approach  to  constructing  and 

235  understanding  microarray  networks  within  a  specific  cell  type  like  the  RPE  of  the 

236  mouse. 

237 

238 

239 

240 

241 

242  Example  of  a  functional  network  in  the  DoD  CDMRP  Retina  Database 

243 

244 

245  To  illustrate  the  features  of  the  new  DoD  CDMRP  Retina  Database,  we  have 

246  chosen  one  specific  gene,  Robo2  (roundabout  homolog  2),  and  will  use  this  gene 

247  to  demonstrate  the  analytical  powers  of  the  database  and  the  bioinformatics  tools 

248  associated  with  GeneNetwork.  Robo2  is  highly  expressed  in  the  retina  with  a 

249  mean  value  of  10.7  across  the  BXD  strain  set.  The  expression  within  individual 

250  strains  varies  from  a  low  of  10.2  to  a  high  of  1 1 .1 .  This  is  a  log2  scale  and 

251  represents  approximately  a  two-fold  difference  in  expression  (Figure  2).  When 

252  we  examine  the  database  for  genes  with  a  similar  pattern  of  expression  across 

253  the  BXD  strain  set,  there  is  a  group  of  genes  that  are  highly  correlated  with  the 

254  expression  pattern  of  Robo2  (Table  3).  One  example  is  the  third  correlate  on  the 

255  list,  Ncam2  (Figure  3)  with  a  value  of  0.926.  Even  the  100th  correlate  on  the  list 

256  {Gitl)  has  a  relatively  high  correlation  (r  =  0.873)  with  Robo2  (See  Supplemental 

257  Table  1). 

258  To  define  the  regions  of  the  genome  modulating  the  expression  of  Robo2,  we 

259  plotted  a  genome  wide  scan  for  Robo2  (Figure  4).  This  plot  defines  regions  of 

260  the  genome  that  correlate  with  the  level  of  Robo2  expression,  a  quantitative  trait 

261  locus  (QTL).  In  this  interval  map  there  is  one  significant  QTL  on  chromosome  16 

262  (notice  the  peak  reaches  the  red  line  on  the  scan,  p  =  0.5)  and  there  are  two 

263  suggestive  peaks  on  Chromosome  1  and  Chromosome  17  (above  the  gray  line). 

264  The  expression  of  Robo2  is  modulated  by  genomic  elements  on  Chromosome 

265  1 6.  There  are  two  types  of  elements  that  could  be  affecting  the  expression  of 

266  Robo2 ;  a  cis-QTL  or  a  gene  with  a  nonsynonymous  SNP.  If  we  examine  the 

267  significant  QTL  on  Chromosome  1 6  (21  -27  Mb),  we  find  there  are  no  significant 

268  cis-QTLs  at  the  gene  level.  With  the  DoD  CDMRP  Retina  Database  it  is  now 

269  possible  to  look  at  the  individual  probes  in  exons  and  introns.  When  we 

270  interrogate  the  DoD  CDMRP  Retina  Exon  Level  Database,  we  find  one  probe 

271  (Affy_1 7329472)  that  lies  within  the  Leprell  gene.  When  we  checked  the 

272  location  of  the  probe  with  the  Verify  function  on  GeneNetwork,  the  probe  lies  in 

273  an  intron  and  may  be  a  non-coding  RNA.  However  when  we  examined  the  RNA- 
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seq  data  from  GeneNetwork  it  appears  that  this  was  detected  in  an  RNA-seq 
analysis  of  the  hippocampus  and  thus  this  may  be  part  of  Leprell  gene  itself. 
Nonetheless,  this  probe  marks  a  candidate  for  modulating  the  expression  of 
Robo2.  The  second  approach  is  to  examine  this  region  for  nonsynonymous 
SNPs.  Using  the  SNP  browser  in  GeneNetwork,  we  looked  at  Chromosome  16 
(21-27  Mb)  and  found  four  known  genes  with  nonsynonymous  SNPs:  Kng2, 

Kngl,  BC106179  and  Maspl.  This  analysis  provides  us  with  five  candidates  for 
modulating  the  expression  of  Robo2. 

To  determine  whether  this  highly  correlated  set  of  genes  in  the  Robo2 
network  have  functional  relationship(s),  we  examined  the  top  500  correlates  of 
Robo2  to  determine  if  there  were  specific  functional  transcript  enrichments  using 
Gene  Set  (WEB-based  GEne  SeT  AbaLysis  Toolkit, 

http://bioinfo.vanderbilt.edu/webgestalt/webgestalt.php).  The  list  of  the  top  500 

correlates  of  F?obo2was  enriched  for  a  number  of  biological  processes  (nervous 
system  development,  synaptic  transmission  and  neuron  differentiations); 
molecular  functions  (enzyme  binding,  PDZ  domain  binding,  inorganic  cation 
transmembrane  transporter,  and  metal  ion  transmembrane  transporter  activity); 
and  cellular  components  (cell  projection  part,  neuron  projection,  intracellular  part, 
and  axon  genes).  This  type  of  analysis  plays  a  critical  role  in  many  genetic 
networks,  defining  the  functional  role  of  the  network.  In  this  specific  case  the 
analysis  demonstrates  that  the  Robo2  network  is  involved  in  axonal  growth  and 
neuronal  development. 

DISCUSSION 

This  paper  announces  the  release  of  two  new  BXD  retina  databases  on 
GeneNetwork.  The  first  is  at  the  gene  level  (DoD  CDMRP  Retina  Affy  MoGene 
2.0  ST  (May15)  RMA  Gene  Level  Database).  The  second  dataset  is  exon  level 
analysis  of  the  same  data  presented  in  the  first  dataset  (DoD  CDMRP  Retina  Affy 
MoGene  2.0  ST  (May15)  RMA  Exon  Level  Database).  In  this  paper  we  attempt  to 
emphasize  some  of  the  special  aspects  of  these  two  datasets  including  an  exon 
level  of  analysis,  the  inclusion  of  microRNAs  and  many  non-coding  RNAs.  To 
illustrate  many  of  these  new  features  we  presented  two  different  approaches  for 
analysis  using  the  datasets. 

The  first  was  an  examination  of  a  cell  signature  within  the  dataset.  Within  the 
DoD  CDMRP  Retina  Database,  there  is  a  pronounced  RPE  signature.  Some 
strains  demonstrate  very  low  levels  of  expression  of  RPE65  while  other  strains 
have  over  16-fold  higher  levels  of  expression.  This  difference  could  not  be  due  to 
differences  in  expression  within  the  RPE,  for  we  know  that  all  RPE  cells  express 
this  gene  at  approximately  the  same  level.  We  believe  that  this  is  due  to 
differences  in  the  time  of  day  the  retinas  were  isolated.  The  retinal  samples  at 
two  different  locations.  At  the  University  of  Tennessee  samples  were  usually 
isolated  starting  at  10:00  AM  and  lights  on  in  the  animal  colony  occurred  at  6:00 
AM.  Thus,  the  retinas  were  isolated  at  least  four  hours  after  lights  on.  At  Emory 
University,  the  retinas  were  isolated  starting  at  9:00  AM  and  lights  on  occurred  at 
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7:00  AM  or  starting  at  2  hours  after  lights  on  in  the  animal  colony.  These 
differences  in  the  time  of  day  the  retinas  were  isolated  may  be  related  to  the 
number  of  RPE  adhering  to  the  retinal  samples  [25]. 


There  are  a  number  of  bioinformatics  tools  available  to  the  vision  research 
community.  These  include  NEI  Bank  project 

(http://neibank.nei.nih.gov/index.shtml ),  which  provides  transcriptome  profiling  of 
the  tissues  of  the  eye  including  mouse  and  human  [26].  The  Cepko  group  has 
provided  the  mouse  retina  serial  analysis  of  gene  expression  (SAGE)  library 
(http://itstqp01.med.harvard.edu/retina)  that  includes  gene  expression  of  the 
embryonic  and  postnatal  retina  [27,  28].  Steve  Dager  and  his  group  have  lists  of 
mapped  loci  and  cloned  genes  associated  with  inherited  retinal  disease  on  the 
RetNet  website  (www.sph.uth.tmc.edu/RetNet/ ).  GENSAT  (Gene  Expression 
Nervous  System  Atlas)  now  has  a  section  devoted  to  the  Retina  Project  [29]  at 
www.qensat.org/retina.isp  .  The  cell  specific  labeling  in  the  retina  for  different 
genes  is  illustrated  using  BAC  transgenic  mice.  The  pattern  of  labeling  in  the 
retina  defines  the  retinal  cells  types  expressing  specific  genes.  This  cellular 
localization  aids  in  defining  localization  of  genetic  networks  in  the  retina.  Finally 
we  have  posted  the  data  from  the  study  of  glaucoma  by  Howell  et  al.[30]  on  the 
GeneNetwork  website  under  the  BXD  eye  database.  These  data  are  very  helpful 
in  understanding  the  role  of  specific  genetic  networks  in  glaucoma  (for  example, 
see  Templeton  et  al.[31  ]). 

In  conclusion,  the  DoD  CDMRP  Retina  Databases  offered  on  GeneNetwork  is 
new  resource  in  an  expanding  variety  of  bioinformatics  tools  available  to  vision 
research  community.  Previously  we  have  offered  several  BXD  microarray 
databases  on  GeneNetwork  included  to  the  vision  science  community:  the 
transcriptome  of  the  whole  eye  (Eye  M430v2  (Sep08)  RMA  Database)  described 
in  detail  by  Geisert  et  al.[18];  a  normal  retina  database  (Normal  Retina  (April 
2010)  Ranklnv  Database)  described  in  detail  by  Freeman  et  al.[19];  and  the 
retina  2  days  after  optic  nerve  crush  (ONC  Retina  (April  2012)  Ranklnv 
Database)  described  in  Templeton  et  al.  [32],  This  new  website  offers  a  unique 
look  at  expression  at  the  exon  level.  In  addition,  there  are  many  non-protein 
coding  transcripts  represented  in  the  dataset.  The  bioinformatics  tools  offered  on 
GeneNetwork  and  these  new  databases  are  a  unique  resource  for  the  vision 
research  community. 
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Figure  1  Expression  of  Rpe65  across  the  BXD  strains  in  the  DoD  CDMRP 
Normal  Retina  Dataset.  The  expression  levels  of  Rpe65  are  shown  for  many  of 
the  BXD  strains  as  mean  expression  and  standard  error  of  the  mean.  The 
individual  strain  identifications  are  shown  along  the  bottom  and  the  scale  is  log2. 
Notice  the  relatively  low  levels  of  Rpe65  in  some  stains  (DBA/2J,  BXD5  BXD12, 
BXD34,  BXD40,  BXD48a,  BXD60,  BXD69,  BXD100,  BXD101  and  BXD102)  and 
8-fold  high  levels  of  expression  other  strains  (BXD1 6,  BXD31 ,  BXD42,  BXD43, 
BXD50,  BXD56,  BXD75  and  BXD85). 
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Figure  2  Expression  of  Robo2  across  the  BXD  strains  in  the  DoD  CDMRP 
Normal  Retina  Dataset.  The  expression  levels  of  Robo2are  shown  for  many  of 
the  BXD  strains  as  mean  expression  and  standard  error  of  the  mean.  The 
individual  strain  identifications  are  shown  along  the  bottom  of  the  plot  and  the 
scale  is  log2.  Notice  the  variability  in  Robo2  expression  across  the  BXD  strains. 
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Figure  3  The  Pearson  correlation  between  Robo2  and  Ncam2.  Ncam2  was  the 
second  highest  correlate  to  Robo2  in  the  DoD  CDMRP.  Notice  that  in  strains 
where  Robo2  is  highly  expressed,  Ncam2  is  also  highly  expressed  and  that  in 
strains  where  Robo2  is  low,  Ncam2  is  also  low. 
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Figure  4.  Genome-wide  Interval  Map  of  Robo2.  This  genome-wide  graph 
displays  the  quantitative  trait  loci  (QTL)  distribution  across  the  DoD  CDMRP 
Normal  Retina  Dataset.  The  X-axis  plots  the  locations  of  the  QTLs  controlling  the 
transcript  expression.  Positions  are  measured  in  mega-bases  from  Chromosome 
1  to  Chromosome  X  (1-2600  Mb).  The  y-axis  plots  Likelihood  Ratio  Statistic 
(LRS).  The  significant  levels  of  individual  QTLs  are  color-coded.  The  red  line 
represents  a  genome-wide  significance  level  (p  >  0.05)  and  the  gray  line  is 
suggestive.  Notice  a  significant  QTL  on  Chromosome  1 6. 


Table  1.  Correlates  of  Rpe65  in  the  DoD  CDMRP  Normal  Retina  Dataset. 
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non-SMC  c  ement  2  homolog 

Chrl  5.  59.205753 

10.909 

0.859 

Plxnal 

plexin  A1 

Chr6:  89  266307 

10.068 

-0.858 

Ubo2p 

ubiq  jilin-con  ugat  ng  enzyme  E2B 

Chrll  51  798648 

11.282 

0.858 

Aacs 

acetoacety  -CoA  synthetase 

Chr5:  125.956184 

9  520 

-0.857 

Acst3 

acyl-CoA  synthelase  ‘amily  member  3 

Ch-8.  125299405 

9  514 

-0.856 

S!c17a7 

soljte  carrier  tarn  ly  1 7,  member  7 

Chr7:  52  419291 

13.143 

-0.856 

Srebf2 

sterol  rogulalu-y  element  binding  laclu-  2 

Chrl  5  81.977696 

11.306 

-0.855 

Epb4.9 

erytnrocyte  protein  band  4.9 

Chrl  4.  71.301070 

9  983 

-0.855 

Vapa 

vusic  e-associaled  membrane  protem.  assocatod  p'otein  A 

Chrl  7:  65.929392 

12.315 

0  855 

Slc12a5 

soljte  came-  lam  iy  12,  member  5 

Ch-2. 164  786302 

12.122 

-0.854 

Card9 

caspase  recruitment  coma  n  family,  memaer  9 

Ch-2:  26  237696 

7.648 

-0.854 

Cfh 

complement  component  lactor  h 

Ch-1:  141  982432 

8.586 

0.853 

Dagla 

diacylglycerol  lipase,  alpha 

Chrl  9.  10.319755 

9.568 

-0.853 

Pde4dip 

phosphodiesterase  4D  interacting  proto  n  (myomegatin) 

Chr3:  97  493751 

9  553 

-0.853 

Ankl 

ankyhn  1 .  eryth-od 

Ch-B:  24.CB5316 

9  543 

-0.853 

Zcchc14 

zinc  linger.  CCFIC  domain  containing  14 

Ch-8:  124.122603 

10.388 

-0.853 

Ss'3 

signal  sequence  receptor,  gamma 

Chr3:  65  186870 

6.036 

0.853 

Mte-lcl 

MTERF  domain  containing  1 

Ch-13:  67.007904 

9547 

0.853 

Apbal 

amylod  beta  (A4)  precursor  protein  bind  ng.  A.  member  1 

Chrl  9.  23  833366 

10.554 

-0.852 

Acot13 

acyl-CoA  tnioeslerase  13 

Ch-13.  24.909817 

10.632 

0.851 

Rbl 

ret  nob  astoma-lixo  1  ;p107) 

Chr2:  156  971629 

9.071 

0.851 

Cpxm2 

carpoxypepldase  X  2  (M14  family) 

Chr7:  139234493 

8  069 

0.851 

Ret 

ret  proto-oncogene 

Ch-6  118.101766 

10.226 

-0.850 

Musya  1 

motiu  spe-m  domain  containing  1 

ChrX:  50.698185 

10.714 

0.850 

Lrp3 

low  dens  ty  1  poprotein  -eceplor-rolalec  proto  n  3 

Chr7.  35.984852 

9  102 

-0.850 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 


Table  2.  MicroRNAs  in  the  Rpe65  signature. 


Rpe65 

Mir98 

Mir449a 

mir301b 

mir28b 

Camk2b 

_ 

_ 

Camk2b 

Atp2b2 

- 

- 

Atp2b2 

- 

Cfl2 

- 

- 

CM2 

- 

Dlcl 

- 

- 

Dlcl 

- 

Eif2c1 

Eif2c1 

- 

Eif2c1 

- 

Nptxl 

Nptxl 

Nptxl 

Nptxl 

- 

Pitpnm2 

- 

— 

Pilpnm2 

- 

Ppp6r1 

- 

- 

Ppp6r1 

— 

Psap 

— 

— 

Psap 

- 

Slc17a7 

Slc17a7 

_ 

Snx2 

_ 

- 

Snx2 

- 

Subl 

- 

— 

Subl 

- 

Tet3 

- 

- 

Tet3 

- 

Zbtb4 

- 

- 

Zbtb4 

- 

Zcchc14 

- 

- 

Zcchc14 

- 

2610507B11Rik 

2610507B11Rik 

2610507B11Rik 

- 

- 

Abr 

Abr 

Abr 

- 

- 

Ahsa2 

Ahsa2 

Ahsa2 

— 

- 

Cacna2C2 

Cacna2d2 

Cacna2d2 

- 

- 

Cntn2 

Cntn2 

Cnln2 

- 

- 

Dcaf7 

Dcaf7 

Dcaf7 

- 

- 

E2f5 

E2f5 

E2f5 

_ 

_ 

FbxolO 

FbxolO 

FbxolO 

- 

- 

Mgat5b 

Mgat5b 

Mgat&b 

_ 

_ 

Ndstl 

Ndstl 

Ndstl 

- 

- 

Nrxn2 

Nrxn2 

Nrxn2 

- 

- 

Pvrll 

Pvrll 

Pvrll 

- 

- 

Ret 

Ret 

Ret 

- 

- 

Slc6a1 

Slc6a1 

Slc6a1 

- 

- 

T(dp2 

Tfdp2 

Tfdp2 

- 

- 

Tnm67 

Trim67 

Tnm67 

- 

- 

Usp31 

Usp31 

Usp31 

- 

- 

Agaol 

Agapl 

- 

- 

- 

Apbal 

Apbal 

- 

- 

- 

Bsn 

Bs'i 

- 

- 

- 

FbxJ14 

FbxJ14 

- 

— 

- 

Insr 

Insr 

_ 

- 

- 

Kcncl 

Kcncl 

- 

- 

- 

Nsmce2 

Nsmce2 

- 

- 

- 

Slc7a14 

Slc7a14 

- 

- 

- 

Srebf2 

Srebf2 

_ 

_ 

_ 

Syt7 

Syt7 

- 

- 

- 

Thbsl 

Thbsl 

- 

- 

- 

Table  3.  Top  20  correlates  of  Robo2. 


Symbol 

Description 

Location  (Chr  Mb) 

Mean  Expr 

8ampler 

Kobo2 

roundabout  nomo  og  2  (Drosophila) 

Ch-16:  73.B92551 

10.72 

1.00 

Casr 

calcs  um 

ChrX:  13.094236 

10.91 

0.94 

Ncam2 

neural  cell  adhesion  molecu  e  2 

Ch-16.  B1.20C942 

9.95 

0.93 

G'ia3 

gl  utamalo  receptor.  ionotropic.  AMPA3 

Ch'X:  3B.754305 

10.25 

0.92 

Lphn3 

latrophlln  3 

Ch-5:  B1  450227 

10.25 

0.92 

Cnks-2 

connocto-  enhancer  o'  Kinase  supp-essor  ol  Has  2 

Ch'X:  154.259368 

10.37 

0.92 

Clstn2 

calsyntonm  2 

Ch-9:  97  344814 

10.00 

0.92 

Gnaq 

guanine  nucootdo  binding  protein 

Ch-19:  16.207321 

10.63 

0.92 

DriajcS 

DnaJ  (Hsp43)  homolog.  suo'amily  C.  momper  6 

Ch'4.  101.169253 

11.97 

0.92 

SIcBal 

solute  earner  lam  ly  8  (socium 

Ch-17.  B1. 772445 

945 

0.91 

Dpys2 

dihydropyrimidinase-  ike  2 

Chr14.  67.421 7C1 

13.25 

0.91 

Plchl 

ptwsohol  paso  C.  eta  1 

Ch'3;  63  500156 

10.96 

0.91 

G'ia1 

glutamate  receptor,  lonot-ooic.  AMPA1  (alpha  1) 

Chrll.  56B24B89 

9.77 

0  91 

Fam165b 

lam  iy  with  sequence  similarity  165,  member  B 

Chr16.  92.301531 

10.04 

-0.91 

N  lu 'ill 

neural  cell  adhesion  molecu  e  1 

Ch-9:  49  310243 

12.95 

0.90 

B4galt6 

UDP-Gal.betaGIcNAc  beta  1  4-galactosyltrans‘eraso 

Ch'18:  2C.B431C0 

10.76 

0.90 

Od2l 

odd  Oc 

ChrX.  40.677132 

8  81 

0.90 

Sortl 

sort  lin  1 

Ch'3:  108.087009 

11.52 

0.90 

Ubqln2 

ubtqullin  2 

ChrX:  149.932775 

11.29 

0.90 

Mtmr9 

myotubularn  related  protein  9 

Chr14:  64.142447 

11.71 

0.90 

Gabbr2 

gamma-am  nobutyrc  aoid  (GABA)  B  receolo-.  2 

Ch-4:  46.675190 

11.04 

0.90 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
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25 

26 

27 

28 

29 

30 

31 

32 
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46  Abstract 

47 

48  Purpose:  The  present  study  was  designed  to  identify  early  changes  associated 

49  with  injury  of  retinal  ganglion  cells. 

50 

51  Methods:  The  normal  retina  database  and  optic  nerve  crush  (ONC)  database  on 

52  GeneNetwork  (www.genenetwork.org)  were  used  to  identify  markers  of  retinal 

53  injury.  One  gene,  Soxll,  was  examined  further  using  two  neuronal  injury 

54  paradigms,  ONC  and  blast  injury  to  the  eye  of  C57BL/6  mice.  The  distribution  of 

55  SOX1 1  was  determined  using  indirect  immunohistochemical  methods  and  the 

56  levels  of  SOX1 1  protein  expression  were  defined  by  semi-quantitative 

57  immunoblot  methods.  In  situ  hybridization  was  performed  using  the  Affymetrix  2- 

58  plex  Quantigene  View  RNA  In  Situ  Hybridization  Tissue  Assay  System. 

59 

60  Results:  SOX11  was  dramatically  upregulated  in  the  retina  following  ONC  and 

61  blast  injury.  The  level  of  Soxll  message  increased  by  approximately  8-fold  2 

62  days  after  ONC.  In  the  normal  retina,  there  was  only  light  immunostaining  for 

63  SOX1 1  in  retinal  ganglion  cells  and  cells  in  the  inner  nuclear  layer.  After  ONC  or 

64  blast  injury,  the  staining  intensity  increased  dramatically  in  both  layers.  In  situ 

65  hybridization  demonstrated  a  similar  distribution  of  message  for  Soxll  in  the 

66  normal  retina  as  well  as  a  profound  increase  in  Soxll  message  within  the 

67  ganglion  cells  following  ONC. 

68 

69  Conclusion:  Taken  together,  these  data  indicate  that  Soxll  may  be  involved  in 

70  the  initial  response  of  the  retina  to  injury,  playing  a  role  in  the  early  attempts  of 

71  the  ganglion  cells  to  survive. 

72 

73 


74 


75  Introduction: 

76 

77  Advances  in  our  ability  to  monitor  molecular  changes  in  neurons  have  led  to  an 

78  increased  understanding  of  the  events  that  transpire  following  neuronal  injury  [1- 

79  5].  After  axonal  damage,  the  initial  response  of  a  central  nervous  system  (CNS) 

80  neuron  may  be  similar  to  that  of  a  neuron  in  the  peripheral  nervous  system  (PNS) 

81  [6].  Ultimately,  neurons  in  the  PNS  will  regenerate  their  axons  and  survive,  while 

82  those  in  the  CNS  do  not  regenerate  their  axons  and  the  cell  bodies  die.  The  initial 

83  response  of  the  CNS  neurons  to  regenerate  and  eventual  failure  of  this 

84  regenerative  response  was  described  by  Ramon  y  Cajal  and  termed  “abortive 

85  regeneration”  [7],  When  examining  the  initial  response  of  neurons  to  injury,  there 

86  appear  to  be  some  common  responses  in  the  CNS  and  PNS.  One  transcription 

87  factor  activated  in  both  the  CNS  and  PNS  after  injury  is  Soxll  [8].  There  is 

88  strong  evidence  that  this  gene  is  part  of  the  transcriptional  network  activated  by 

89  injury  and  involved  in  axonal  regeneration  in  the  PNS  [9,  10]. 

90 

91  Soxll  is  a  member  of  the  SRY-related  box  group  C  ( SoxC )  gene  family  of 

92  transcription  factors  [11,  12],  SOX11  along  with  SOX4  play  a  critical  role  in  the 

93  normal  development  of  neurons  and  specifically  retinal  ganglion  cells  [13-15]. 

94  SOX1 1  is  expressed  in  retinal  progenitor  cells  as  part  of  the  process  leading 

95  progenitor  cells  to  become  neuroblasts  [16,  17],  In  knock-downs  of  either  Soxll 

96  or  Sox4,  there  is  a  moderate  reduction  in  retinal  ganglion  cell  number;  however 

97  when  both  Soxll  and  Sox4  are  knocked-down,  there  is  a  complete  loss  of 

98  ganglion  cell  development  [13].  During  eye  development,  Soxll  is  also  required 

99  to  maintain  proper  levels  of  hedgehog  signaling,  and  mutations  have  been 

100  associated  with  coloboma  due  to  improper  optic  fissure  closure  [18,  19]. 

101  Furthermore,  SOX11  is  critical  for  axonal  growth,  driving  the  expression  of  axon 

102  growth-related  proteins  such  as  class  III  beta  tubulin  and  MAP2  [20].  SOX1 1  also 

103  plays  a  similar  role  in  adult  neurogenesis.  High  levels  of  SOX1 1  are  found  in  the 

104  cells  within  the  subventricular  zone,  the  rostral  migratory  stream  and  within  the 


105  neuroprogenitor  zone  of  the  dentate  gyrus  [15,  16,  21].  These  studies  underline 

106  the  importance  of  SOX1 1  in  terminal  differentiation  of  progenitor  cells  to  neurons 

107  and  axon  extension. 

108 

109  In  addition  to  functioning  in  neuronal  differentiation,  SOX11  has  a  prominent  role 

110  in  the  response  of  neurons  to  injury.  After  peripheral  nerve  injury,  SOX11  is 

in  immediately  upregulated  in  the  neuronal  cell  bodies  as  the  axon  is  regenerating 

112  [9,  22],  Decreasing  levels  of  SOX11  in  the  neuronal  cell  body  results  in  slower 

113  axonal  regeneration  of  peripheral  nerves  [9j.  Similar  results  are  observed  in 

114  tissue  culture.  When  Soxll  is  knocked  down  in  cultured  peripheral  neurons, 

115  there  is  also  a  reduction  in  neurite  growth  and  an  increase  in  apoptosis  [23]. 

116  Conversely,  overexpressing  Soxll  in  cultured  dorsal  root  ganglion  cells 

117  produces  an  increase  in  neurite  growth,  and  in  vivo  overexpression  of  Soxll 

118  accelerates  the  growth  of  regenerating  axons  [10].  One  intriguing  anatomical 

119  experimental  model  is  the  dorsal  root  ganglion,  where  the  central  projection  of 

120  the  dorsal  root  ganglion  enters  the  spinal  cord  (CNS)  and  the  peripheral 

121  projection  extends  out  into  a  peripheral  nerve  that  is  myelinated  by  Schwann 

122  cells.  When  the  central  rootlet  is  severed,  there  is  a  modest  (51%)  increase  in 

123  Soxll  expression  in  the  ganglion  even  when  the  central  portion  will  not 

124  regenerate  back  into  the  spinal  cord.  However,  when  the  peripheral  root  is 

125  damaged,  a  relatively  massive  (1004%)  increase  in  Soxll  is  seen  as  the  axons 

126  regenerate  down  the  peripheral  nerve. 

127  In  the  present  study,  we  examine  the  role  of  SOX11  in  the  retina  following 

128  injuries  to  the  axons  of  the  optic  nerve.  We  propose  that  the  upregulation  of 

129  SOX11  after  injury  is  an  attempt  of  neurons  to  regenerate,  but  ultimately  results 

130  in  abortive  regeneration  and  cell  death. 

131 


132  Materials  and  Methods 

133 

134  Optic  Nerve  Crush.  The  Optic  Nerve  Crush  procedure  was  performed  as 

135  previously  described  in  Templeton  et  al  [24],  Briefly,  C57BL/6  (n  =  5)  mice  were 

136  deeply  anesthetized  with  13  mg/kg  Rompum  and  87  mg/kg  Ketalar  for  the 

137  surgery.  An  incision  was  made  into  the  lateral  aspect  of  the  conjunctiva,  and  the 

138  eye  was  rotated  nasally  to  expose  the  optic  nerve.  The  optic  nerve  was  then 

139  grasped  for  10  seconds  with  Dumont  cross-clamp  #7  forceps  (Roboz,  cat. 

140  #RS=5027,  Gaithersburg,  MD),  using  only  the  spring  action  of  the  instrument  to 

141  crush  the  nerve.  The  instrument  was  then  removed  and  the  eye  was  allowed  to 

142  rotate  back  into  place.  The  animals  were  allowed  to  recuperate  from  surgery  on  a 

143  water-heated  warming  pad.  The  Institutional  Animal  Care  and  Use  Committee 

144  approved  all  procedures  for  the  surgery  and  handling  of  mice  (at  the  University  of 

145  Tennessee  Health  Science  Center  and  Emory  University). 

146 

147  Blast  Injury  to  the  Eye.  Blast  injuries  to  the  eye  were  produced  as  previously 

148  reported  in  Hines-Beard  et  al  [25].  Briefly,  the  mice  were  deeply  anesthetized  and 

149  secured  in  a  PVC  pipe.  A  50psi  overpressure  wave  was  delivered  selectively  to 

150  the  right  eye  of  C57BL/6  mice  (n  =  5).  The  output  pressure  of  the  blast  apparatus 

151  was  measured  immediately  before  and  after  the  procedure,  ensuring  an  accurate 

152  and  precise  injury.  The  animals  were  sacrificed  5  days  after  the  blast  injury.  The 

153  Institutional  Animal  Care  and  Use  Committee  approved  all  procedures  for  the 

154  surgery  and  handling  of  mice  (at  the  University  of  Tennessee  Health  Science 

155  Center  and  Emory  University). 

156 

157  Immunoblot  Analysis.  Protein  samples  from  12  retinas  were  examined  to 

158  determine  the  relative  levels  of  SOX11.  Animals  were  anesthetized  with  a 

159  mixture  of  xylazine  (13  mg/kg,  Rompum)  and  ketamine  (87  mg/kg,  Ketalar)  and 

160  following  a  cervical  dislocation  the  retinas  were  removed.  The  dissected  retinas 

161  were  placed  in  IX  reducing  sample  buffer  (2%  SDS  and  10%  glycerol  in  0.05  M 

162  Tris-HCI  buffer,  pH  6.8).  Equal  amounts  of  protein  were  loaded  onto  each  lane  of 


163  an  SDS  PAGE  gel.  The  balance  of  the  proteins  was  determined  by  coomassie 

164  blue  staining  and  subsequent  quantification  of  total  protein  in  each  lane.  To 

165  analyze  the  concentrations  of  SOX11  in  the  retina,  the  proteins  were  transferred 

166  to  PVD  membranes.  The  blots  were  blocked  with  2%  non-fat  dry  milk  in 

167  phosphate  buffer  (pH  7.4)  and  probed  overnight  with  the  rabbit  anti-  SOX11 

168  primary  antibody  (Santa  Cruz  Biotechnology,  Inc.  California).  We  then  rinsed  the 

169  blots  and  probed  them  with  the  HRP-labeled  donkey  anti-rabbit  secondary 

170  antibody  (Jackson  ImmunoResearch  Laboratories,  Inc.,  West  Grove, 

171  Pennsylvania).  The  blots  were  rinsed  with  0.5  M  Tris  buffer  (pH  7.4)  and  reacted 

172  with  0.05%  DAB  and  hydrogen  peroxide.  The  blots  were  washed  and  visualized 

173  using  an  ECL  detection  kit  (Thermo  Fisher  Scientific  Inc.  Rockford,  Illinois)  and 

174  Kodak  4000  MM  image  station.  To  define  the  amount  of  protein  loaded  onto  each 

175  lane  the  blot  was  stripped  and  restained  for  actin  using  a  mouse  anti-B-Actin 

176  antibody  (A5441,  Sigma  Aldrich,  St  Louis  MO)  followed  by  an  HRP-labeled 

177  Donkey  anti-mouse  secondary  antibody  (Jackson  ImmunoResearch 

178  Laboratories,  Inc.,  West  Grove,  Pennsylvania).  The  blots  were  then  rescanned. 

179 

180  Immunohistochemistry.  For  the  immunohistochemical  analysis  of  the 

181  distribution  of  SOX11,  ten  C57BL/6  mice  (2  control  mice  and  2  mice  of  each 

182  group  at  both  2  and  5  days  postoperatively)  were  used.  The  mice  were  deeply 

183  anesthetized  with  2,2,2-tribromoethanol  and  then  perfused  transcardially  with 

184  saline  followed  by  4%  paraformaldehyde  using  a  peristaltic  pump  (Cole  Parmer 

185  Instruments,  Chicago,  IL).  After  an  hour  of  post-fixation  in  paraformaldehyde  at 

186  room  temperature,  retinas  were  rinsed  in  0.1  M  phosphate-buffered  saline  (PBS), 

187  embedded  in  4%  low-melt  agarose,  and  cut  on  a  vibratome  into  50  pm  thick 

188  slices.  Sections  were  then  blocked  in  4%  bovine  serum  albumin  (Sigma-Aldrich, 

189  St.  Louis,  MO)  in  PBS  with  0.1%  Triton  X-100  (Sigma)  on  a  rotating  platform  for  1 

190  hour,  washed  3  times  and  incubated  with  primary  antibody  at  4°C  overnight.  The 

191  next  day,  sections  were  rinsed  again  and  then  incubated  in  Alexa  Fluor- 

192  conjugated  secondary  antibody  (Life  Technologies).  After  that,  sections  were 

193  rinsed  in  PBS,  mounted  with  Fluoromount  (Southern  Biotech,  Birmingham,  AL) 


194  and  coverslipped.  The  following  primary  antibodies  and  dilutions  were  used:  anti- 

195  SOX11  (US  Biological  S5364-40C,  1:500),  anti-TUJI  (Tujl,  a  gift  from  Anthony 

196  Franfurter  [26]),  and  anti-PKC  alpha  (abeam  ab32376,  1 :250).  The  following  anti- 

197  SOX11  antibodies  were  tried  for  Immunohistochemistry  but  did  not  work:  Santa 

198  Cruz  sc-20096,  abeam  ab170916,  biorbyt  101332,  and  LifeSpan  Biosciences  LS- 

199  Cl  0306. 

200 

201  Confocal  Microscopy.  The  sections  were  scanned  and  10  pm  thick  z-stacks 

202  were  acquired  on  a  Eclipse  Ti  confocal  microscope  (Nikon  Inc.,  Melville,  NY) 

203  equipped  with  Nikon  Cl  software.  Scans  were  taken  using  a  40x  Plan  Fluor 

204  40x/1 .30/0.20  Oil  objective.  Laser  power  and  gain  were  kept  constant  for  all 

205  pictures.  Images  were  loaded  into  Fiji  [27]  and  average  intensity  projections  were 

206  created.  The  final  composition  of  the  images  was  made  using  Photoshop  CS  6 

207  (Adobe  Systems  Inc.,  San  Jose,  CA). 

208 

209 

210  In  Situ  Hybridization.  In  situ  hybridization  was  performed  using  the  2-plex 

211  Quantigene  View  RNA  ISH  Tissue  Assay  kit  (Affymetrix,  Inc.,  Santa  Clara, 

212  California  ).  Assays  were  performed  as  per  the  manufacturer’s  instructions  with 

213  stock  solutions.  Eyes  were  isolated  from  C57BL/6  mice  (n  =  6),  both  control  and 

214  two  days  after  optic  nerve  crush,  and  drop-fixed  in  4%  Paraformaldehyde  for  24 

215  hours.  Mid-way  through  the  24-hour  fixation,  a  26-gauge  needle  was  used  to 

216  create  a  hole  to  the  vitreous  humor,  assisting  with  the  fixation.  Immediately  after 

217  the  24-hour  period,  the  eyes  were  serially  dehydrated  in  ethanol  and  embedded 

218  in  paraffin.  Blocks  were  sectioned  on  an  American  Optical  series  1000  microtome 

219  (American  Optical  Co.  Buffalo,  New  York)  to  a  thickness  of  5  pm  and  mounted  on 

220  Surgipath  X-tra  micro  slides  (Leica  Biosystems  Richmond  Inc.,  Richmond, 

221  Illinois).  Before  beginning  the  hybridization  protocol,  the  slides  were  baked  at  60° 

222  C  for  30  minutes  to  increase  tissue  adhesion.  As  per  the  manufacturer’s  protocol, 

223  the  paraffin  was  removed  from  the  slides  with  xylene  before  being  boiled  in  a 

224  pretreatment  solution  (Affymetrix)  for  10  minutes  and  incubated  with  Protein 


225  Kinase  K  [Affymetrix  Santa  Clara,  CA]  at  40°  C  for  10  minutes.  Custom  probes  for 

226  Soxll  and  Chrna6  [Affymetrix  Santa  Clara,  CA]  were  then  hybridized  to  the 

227  tissue.  Signal  amplification  was  accomplished  by  hybridizing  Type  1  and  Type  2 

228  specific  pre-amp  oligonucleotides,  amp  oligonucleotides  and  label 

229  oligonucleotides  sequentially,  achieving  a  400-fold  signal  amplification  from  each 

230  mRNA  molecule.  Sections  were  viewed  with  an  Olympus  BX51  microscope 

231  (Olympus  America  Inc.,  Melville,  New  York). 

232 


233  Results 

234  We  first  compared  Soxll  mRNA  levels  in  normal  mice  to  Soxll  levels  following 

235  ONC  using  the  bioinformatic  tools  on  GeneNetwork  (genenetwork.org)  and  two 

236  preexisting  databases  generated  by  us:  Normal  HEI  Retina  (April  2010)  and  ONC 

237  HEI  Retina  (April.  2012)  [4,  28].  Soxll  was  one  of  the  genes  with  the  largest 

238  change  in  expression  two  days  after  optic  nerve  crush.  In  the  normal  retinal 

239  dataset,  the  mean  expression  for  Soxll  (detected  by  lllumina  probe 

240  I LMN_1 235647)  across  the  BXD  Rl  strains  was  8.4  on  a  2  Z  +  8  Log2  scale  (this 

241  is  just  above  the  mean  detection  level  of  mRNA  on  the  array,  which  is  set  to  8). 

242  In  the  C57BL/6  parental  strain,  the  expression  level  in  the  normal  retina  was  8.59 

243  and  for  the  DBA/2J  strain  the  mean  expression  level  was  8.54.  Two  days  after 

244  ONC,  there  was  a  dramatic  increase  in  the  level  of  Soxll  expression  (Figure  1), 

245  with  the  mean  expression  across  the  BXD  strains  being  11.03,  which 

246  corresponded  to  an  approximately  8-fold  increase.  The  same  increase  was 

247  observed  in  individual  strains.  The  C57BL/6  strain  had  an  expression  level  of 

248  11.33  after  ONC  and  the  expression  in  the  DBA/2J  strain  increased  to  11.44. 

249  These  data  indicate  that  Soxll  is  dramatically  upregulated  after  a  specific  injury 

250  to  the  ganglion  cell  axons  within  the  optic  nerve. 

251 

252  In  situ  hybridization  was  used  to  identify  the  cells  expressing  Soxl  1  after  injury  to 

253  the  retina.  For  the  in  situ  hybridization,  we  examined  retinas  two  days  after  optic 

254  nerve  crush  (Fig.  2B)  and  compared  it  to  uninjured  control  retinas  (Fig.  2A). 

255  Using  the  Affymetrix  2-plex  Quantigene  View  RNA  ISH  Tissue  Assay  kit,  we 

256  labeled  cells  expressing  the  retinal  ganglion  cell  marker  Chrna6  blue  and  Soxll 

257  red  (Fig.  2)  [3j.  In  the  control  retina,  many  of  the  cells  in  the  ganglion  cell  layer 

258  were  heavily  labeled  for  Chrna6,  and  a  few  of  the  cells  expressed  low  levels  of 

259  Soxll.  Two  days  after  optic  nerve  crush  there  was  a  dramatic  decrease  in  the 

260  labeling  for  Chrna6  and  a  substantial  increase  in  the  amount  of  labeling  for 

261  Soxll.  This  was  similar  to  the  changes  in  message  levels  observed  in  our 

262  microarray  databases  on  GeneNetwork.org.  In  the  normal  retina  database, 

263  Chrna6  (probe  ILMN_2732438)  was  expressed  at  relatively  high  levels  (mean 


264  value  of  11.04)  and  two  days  following  optic  nerve  crush,  the  expression 

265  decreased  two-fold  (mean  value  of  9.96).  These  data  from  the  in  situ 

266  hybridization  mirror  our  microarray  expression  analysis  results. 

267 

268  Indirect  immunohistochemistry  was  used  to  define  the  distribution  of  SOX1 1 

269  protein  in  the  retina.  To  identify  the  cell  types  expressing  SOX11  after  retinal 

270  injury,  we  stained  the  retina  for  SOX11  and  for  class  III  beta  tubulin,  a  known 

271  ganglion  cell  marker  [26].  In  these  double-labeled  sections,  staining  for  SOX11 

272  was  colocalized  with  class  III  beta  tubulin  in  retinal  ganglion  cells  (Fig.  3), 

273  indicating  that  they  were  also  positive  for  SOX1 1 .  Both  markers  labeled  the  same 

274  population  of  cells  in  the  ganglion  cell  layer  and  the  labeling  extended  out  into 

275  dendritic  processes.  The  staining  pattern  also  indicated  that  SOX11  was  present 

276  in  the  cytoplasm  and  the  nucleus.  Even  though  nuclear  labeling  was  light  in 

277  comparison,  there  were  many  cases  where  a  negative  spot  (presumably  the 

278  nucleolus)  could  be  observed.  These  data  demonstrate  that  SOX1 1  is  present  in 

279  injured  retinal  ganglion  cells. 

280 

281  Using  immunohistochemistry,  we  also  examined  the  change  in  expression  of 

282  SOX1 1  following  optic  nerve  crush  (Fig  4B  and  4C)  and  a  blast  injury  to  the  eye 

283  (Fig.  4E  and  4F).  In  both  cases,  there  was  a  noticeable  increase  in  ganglion  cell 

284  staining  relative  to  that  observed  in  the  uninjured  control  section  (Fig.  4A).  The 

285  staining  was  primarily  cytosolic  and  extended  out  into  the  dendritic  processes  of 

286  the  cells.  In  the  normal  retina,  we  could  identify  light  labeling  of  cells  within  the 

287  inner  nuclear  layer.  The  intensity  of  labeling  within  the  retinal  ganglion  cells 

288  increased  dramatically  in  the  injured  retinas  (Fig.  4B,  4C,  4E  and  4F),  and  there 

289  was  also  an  increase  in  staining  of  cells  within  the  inner  nuclear  layer,  which 

290  colocalized  partly  to  bipolar  cells  (Fig.  5). 

291 

292  The  upregulation  of  SOX1 1  following  ONC  was  confirmed  with  Immunoblotting. 

293  The  intensity  of  the  SOX11  band  from  retinas  after  nerve  crush  (Fig.  6,  lanes  A 

294  and  C)  is  higher  than  in  the  control  retina  samples  (Fig.  6  lanes  B  and  D).  Thus, 


295  there  is  a  dramatic  increase  in  the  expression  of  both  Soxl  1  mRNA  and  SOX1 1 

296  protein  following  retinal  injury. 

297 


298  Discussion 

299  In  the  present  study,  we  demonstrate  that  SOX1 1  is  upregulated  following  retinal 

300  injury.  For  retinal  ganglion  cells,  there  are  many  known  markers,  such  as: 

301  Chrna6,  Pou4f1,  Tubb3,  Thyl,  and  Sncg  [3,  13,  29-33].  Using  these  marker 

302  genes,  we  performed  a  meta-analysis  of  their  expression  in  injured  and  uninjured 

303  retinas  (data  not  shown).  Our  injury  paradigms  included  controlled  optic  nerve 

304  crush  and  selective  blast  injury  to  the  eye,  procedures  that  are  known  to  induce 

305  retinal  ganglion  cell  apoptosis  [24,  25].  As  expected,  the  expression  of  these 

306  retinal  ganglion  cell  markers  decreases  as  the  ganglion  cells  die.  Our  data 

307  indicate  an  overall  reduction  in  expression  of  these  markers,  both  2  and  5  days 

308  after  retinal  insult  [34],  Specifically,  a  3-fold  reduction  in  Sncg  expression  is 

309  observed,  Chrna6  and  Pou4f1  are  down  regulated  2-fold,  and,  while  not  as 

310  robust,  Thyl  and  Tubb3  also  show  a  decrease  in  expression.  Collectively,  these 

311  data  demonstrate  that  known  retinal  ganglion  cell  markers  decrease  as  the  optic 

312  nerve  degenerates. 

313 

314  Unlike  these  markers  for  retinal  ganglion  cells,  Soxll  levels  increase 

315  substantially  following  injury  to  the  retina  and  optic  nerve.  Two  days  after  ONC, 

316  there  was  an  almost  8-fold  upregulation  of  Soxll.  We  wondered  if  this  increase 

317  in  expression  was  specific  to  the  severity  of  retinal  injury,  and  queried  a  publicly 

318  available  microarray  dataset  of  glaucomatous  mice  presented  by  Howell  and 

319  colleagues  as  an  independent  test  on  the  role  of  Soxll  in  the  response  of  the 

320  retina  to  injury  [1],  In  this  DBA/2J  mouse  model  of  pigment  dispersion  glaucoma, 

321  there  were  relatively  low  levels  of  Soxll  in  the  control  mice  and  in  mice  that  did 

322  not  have  detectable  levels  of  ganglion  cell  loss  (Fig.  7).  However,  in  mice  with 

323  moderate  ganglion  cell  loss,  levels  of  Soxll  were  approximately  four-fold  higher 

324  than  in  control  animals.  In  animals  with  severe  glaucoma,  the  level  of  Soxll 

325  decreased  to  near  control  levels.  These  data  suggest  that  the  levels  of  Soxll 

326  decrease  as  the  ganglion  cells  die.  This  transient  upregulation  of  Soxll  in 

327  moderate  cases  demonstrates  the  exquisite  ability  of  Soxll  to  mark  injured, 


328  potentially  dying  neurons.  As  such,  our  data  indicate  that  Soxll  is  a  novel 

329  marker  for  injured  neurons,  specifically  retinal  ganglion  cells. 

330 

331  What  is  the  potential  role  of  SOX1 1  following  injury  to  the  axons  of  the  retinal 

332  ganglion  cell?  One  hint  comes  from  studies  on  the  response  of  dorsal  root 

333  ganglion  neurons  to  peripheral  nerve  injury.  Soxll  is  dramatically  upregulated  in 

334  the  dorsal  root  ganglion  following  injury  to  the  peripheral  nerve  and  plays  a 

335  pivotal  role  in  axonal  regeneration  [9],  Three  days  after  the  transection  of  the 

336  sciatic  nerve  in  the  rat,  there  is  a  1004%  increase  in  Soxll  in  the  dorsal  root 

337  ganglion  and  these  levels  remain  elevated  for  at  least  the  next  4  days  [9].  By  4 

338  weeks  after  transection,  Soxll  levels  have  returned  to  baseline.  This 

339  upregulation  and  sustained  expression  of  Soxll  is  critical  to  the  survival  of  the 

340  dorsal  root  ganglion  neurons  and  the  regeneration  of  peripheral  axons  along  the 

341  injured  nerve.  When  the  levels  of  Soxll  were  knocked  down  by  delivery  of 

342  siRNA,  there  was  an  increase  in  apoptosis  in  cultured  neurons  as  well  as  a 

343  decrease  in  neurite  outgrowth  in  culture  [23].  A  similar  decrease  in  axonal 

344  regeneration  occured  in  vivo  after  Soxll  knockdown  in  the  dorsal  root  ganglion 

345  [9].  Thus,  the  upregulation  of  Soxll  appears  to  be  necessary  for  the  normal 

346  regeneration  of  axons  in  the  peripheral  nerve.  Interestingly,  Soxll  was  also 

347  upregulated  when  the  central  projection  of  the  dorsal  root  was  severed;  however, 

348  the  degree  of  upregulation  was  lower,  only  51%  [9,  23].  Just  like  the  neurons  in 

349  the  retina,  the  dorsal  root  ganglia  did  not  regenerate  their  central  projections  into 

350  the  adjacent  spinal  cord  [35,  36].  This  begs  the  question,  why  is  Soxll 

351  upregulated  in  neurons  of  the  retina  and  the  dorsal  root  ganglion,  when  the 

352  axons  will  not  successfully  regenerate?  The  most  parsimonious  explanation  is 

353  that  these  neurons  are  attempting  to  survive  using  the  same  program  that  is 

354  successful  in  peripheral  nerve  injury;  however,  other  influences  derail  the 

355  regenerative  program  causing  the  cell  to  abort  the  regenerative  process  and,  in 

356  many  cases,  cause  subsequent  neuronal  death  [37-40].  Recent  evidence 

357  indicates  that  these  inhibitory  influences  can  at  least  be  partially  averted  by 

358  stimulating  the  ganglion  cells  to  regrow  their  axons  [41-47],  Given  that 


359  upregulation  of  Soxl  1  is  among  the  initial  responses  of  the  neuron  to  injury  and 

360  that  this  response  is  intact  in  retinal  ganglion  cells,  it  may  be  possible  to  identify 

361  where  the  transcriptional  cascade  leading  to  axon  regeneration  and  cell  survival 

362  in  the  CNS  differs  from  that  in  peripheral  nerve  injury. 

363 
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Figure  1.  Soxll  Expression  in  the  Normal  Retina  and  After  Optic  Nerve 
Crush.  The  ordinate  represents  mRNA  levels  from  microarrays,  expressed  in  a 
log2  scale  with  the  mean  set  to  8.  The  mice  used  to  generate  these  data  were 
C57BL/6J,  DBA/2J,  their  respective  FI  crosses,  BALB/cByJ,  and  members  of  the 
BXD  recombinant  inbred  strain  line.  The  top  graph  denotes  the  expression  levels 
of  Soxll  under  normal  conditions,  while  the  bottom  graph  denotes  the 
expression  of  Soxll  two  days  after  the  optic  nerve  crush.  Note  the  dramatic 
increase  in  expression  following  injury,  from  a  normal  mean  expression  value  of 
8.4  to  a  mean  expression  value  of  1 1 .0  after  nerve  crush. 
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Figure  2.  In  Situ  Hybridization  for  Soxll  and  Chrna6.  The  distribution  of 
Soxll  and  Chrna6  mRNA  was  defined  using  two-coior  in  situ  hybridization  for  a 
section  of  control  retina  (A)  and  a  section  of  a  retina  2  days  after  optic  nerve 
crush  (B).  The  probe  for  Chrna6  is  shown  in  blue  and  that  for  Soxll  in  red.  Note 
the  extensive  staining  for  Chrna6  within  the  ganglion  cell  layer  in  the  control 
section  (Arrow  in  A).  Two  days  after  optic  nerve  crush  there  is  already  a 
substantial  decrease  in  Chrna6  expression.  Very  low  levels  of  Soxll  are  seen  in 
the  control  retina.  There  is  a  dramatic  increase  in  Soxl  1  staining  after  nerve 
crush  within  the  retinal  ganglion  cells  (Arrow  in  B).  A  and  B  are  taken  at  the 
same  magnification  and  the  scale  bar  in  B  =  25  pm. 


531 

532  Figure  3.  S0X11  Expression  in  Retinal  Ganglion  Cells.  A  section  of  retina, 

533  two  days  after  optic  nerve  crush,  was  double  stained  for  SOX1 1  (A)  and  class  III 

534  beta  tubulin  (B),  using  double  label  indirect  immunohistochemistry.  Class  III  beta 

535  tubulin  is  a  marker  for  retinal  ganglion  cells.  In  the  merged  image  (C),  retinal 

536  ganglion  cells  are  double  labeled  (arrowhead)  with  SOX11  and  Class  III  beta 

537  tubulin.  Scale  bar  in  C  =  20  pm. 
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Figure  4.  Upregulation  ofSOXII  After  Retinal  Injury. 

Retinal  sections  were  stained  with  an  antibody  directed  against  SOX1 1  in  the 
following  conditions:  normal  retina  (A),  retina  two  (B)  and  five  (C)  days  after  optic 
nerve  crush,  and  retina  two  (E)  and  five  (F)  days  after  blast  injury.  The  negative 
control  without  primary  antibody  is  shown  in  (D).  After  injury  to  the  retina,  there  is 
a  dramatic  increase  in  the  intensity  of  staining  within  the  cell  bodies  of  retinal 
ganglion  cells,  and  some  cells  in  the  INL.  The  layers  of  the  retina  are  labeled: 
ONL  (outer  nuclear  layer),  INL  (inner  nuclear  layer)  and  GCL  (ganglion  cell 
layer).  All  photomicrographs  are  presented  at  the  same  magnification  and  the 
scale  bar  in  F  =  20  pm. 
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Figure  5.  Upregulation  of  SOX11  in  the  Inner  Nuclear  Layer  Following 
Injury. 

Retinas  were  stained  with  antibodies  directed  against  PKC  alpha  to  label  bipolar 
cells  (A)  and  Soxll  (B)  two  days  after  optic  nerve  crush.  In  the  merged  image 
(C),  about  half  of  the  bipolar  cells  are  double-labeled,  and  they  are  mostly 
localized  to  the  inner  part  of  the  INL  (arrowheads). 

The  layers  of  the  retina  are  labeled:  ONL  (outer  nuclear  layer),  INL  (inner  nuclear 
layer)  and  GCL  (ganglion  cell  layer).  All  photomicrographs  are  presented  at  the 
same  magnification  and  the  scale  bar  in  C  =  20  pm. 
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566  Figure  6.  Upregulation  of  SOX11  protein  after  optic  nerve  crush. 

567  Immunoblots  of  protein  samples  from  normal  retina  (B  and  D)  and  from  retinas  2 

568  days  after  optic  nerve  crush  (A  and  C)  were  probed  with  antibodies  directed 

569  against  SOX1 1.  There  was  increased  staining  of  the  bands  from  the  optic  nerve 

570  crush  samples  relative  to  the  normal  retinal  samples.  Loading  control  in  E-H: 

571  beta-actin. 
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Figure  7.  Soxll  mRNA  Levels  in  Pigmentary  Glaucoma.  Soxll  levels  in  a 
meta-analysis  of  a  DBA/2J  mouse  pigmentary  glaucoma  dataset  from  Howell  et 
al.  (2011).  The  ordinate  represents  the  Soxll  mRNA  level  from  microarrays, 
expressed  in  log2,  scaled  with  the  mean  set  to  8.  The  mice  were  classified  as 
wild  type  Gpnmb+,+  controls  (Gpnmb+  Controls),  no  detectable  glaucoma  1  (NOE 
1),  no  detectable  glaucoma  2  (NOE2),  moderate  glaucoma  (Moderate),  and 
severe  glaucoma  (Severe).  Other  than  the  wild  type  group,  all  the  groups  were 
Gpnmb~'~.  Notice  that  there  is  a  dramatic  increase  in  Soxll  expression  during 
early  phases  of  glaucoma  in  this  mouse  model.  As  the  ganglion  cells  die  in 
severe  glaucoma,  the  expression  levels  of  Soxll  appear  to  decrease.  These 
data  suggest  that  Soxll  is  an  early  marker  for  retinal  injury.  This  general 
expression  pattern  was  observed  for  the  following  5  of  the  6  Affymetrix  probe 
sets  targeting  Soxll :  1429051_s_at  (illustrated  above),  1453002_at, 

1 429372_at,  1453125_at  and  1431225_at. 


